CP Violation in a Supersymmetric S'O(IO) x U{2)f Model 
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A model based on SUSY SO(IO) combined with (7(2) family symmetry constructed recently by 
the authors is generalized to include phases in the mass matrices leading to CP violation. In contrast 
with the commonly used effective operator approach, 126-dimensional Higgs fields are utilized to 
construct the Yukawa sector. _R-parity symmetry is thus preserved at low energies. The symmetric 
mass textures arising from the left-right symmetry breaking chain of 50(10) give rise to very good 
predictions for quark and lepton masses and mixings. The prediction for sin 2/3 agrees with the 
average of current bounds from BaBar and Belle. In the neutrino sector, our predictions are in good 
agreement with results from atmospheric neutrino experiments. Our model accommodates both the 
LOW and QVO solutions to the solar neutrino anomaly; the matrix element for neutrinoless double 
beta decay is highly suppressed. The leptonic analog of the Jarlskog invariant, Jcp, is predicted to 
be of 0(10"^). 

PACS numbers: 12.15Ff,12.10Kt,14.60Pq 



5*0(10) has long been thought to be an attractive can- 
didate for a grand unified theory (GUT) for a number of 
reasons: First of all, it unifies all the 15 known fermions 
with the right-handed neutrino for each family into one 
16-dimensional spinor representation. The seesaw mech- 
anism then arises very naturally, and the non-zero neu- 
trino masses can thus be explained. Since a complete 
quark-lepton symmetry is achieved, it has the promise 
for explaining the pattern of fermion masses and mixing. 
Because B—L contained in 50(10) is broken in symmetry 
breaking chain to SM, it also has the promise for baryoge- 
nesis. Recent atmospheric neutrino oscillation data from 
Super-Kamiokande indicates non-zero neutrino masses. 
This in turn gives very strong support to the viability 
of 50(10) as a GUT group. Models based on 50(10) 
combined with discrete or continuous family symmetry 
have been constructed to understand the flavor prob- 
lem ^ . Most of the models utilize "lopsided" mass tex- 
tures which usually require more parameters and there- 
fore are less constrained. Furthermore, the right-handed 
neutrino Majorana mass operators in most of these mod- 
els are made out of 16h x 16h which breaks the i?-parity 
at a very high scale. We have recently constructed a re- 
alistic model based on supersymmetric 50(10) combined 
with U (2) family symmetry [Q (referred to as "CM" here- 
after) which successfully predicts the low energy fermion 
masses and mixings. Since we utilize symmetric mass 
textures and 126-dimensional Higgs representations for 
the right-handed neutrino Majorana mass operator, our 
model is more constrained in addition to having i?-parity 
conserved Q . The aim of this paper is to generalize this 
model to include phases in the mass matrices which lead 
to CP violation. We first summarize our model followed 
by analytic analyses of the complex mass textures. And 
then the implications of our model for neutrino mixing 



and CP violation are presented. 

The Model: The details of our model based on 
50(10) X U{2)f are contained in CM. The following is 
an outline of its salient features. In order to specify the 
superpotential uniquely, we invoke Z2 x Z2 x Z2 discrete 
symmetry. The matter fields are 



V'a~(16,2)- 



(a = 1,2), ^3^(16,1) 



+++ 



where a = 1, 2 and the subscripts refer to family indices; 
the superscripts +/— refer to (^2)'^ charges. The Higgs 
fields which break 50(10) and give rise to mass matrices 
upon acquiring VEV's are 



(10,1) 
(126,1) 



T2 



T4 



(1) 



Higgs representations 10 and 126 give rise to Yukawa cou- 
plings to the matter fields which are symmetric under the 
interchange of family indices. The left-right symmetry 
breaking chain of 50(10) is 



50(10) — * 5C/(4) X SU{2)l x SU{2)r 

5[/(3) X SU{2)l X SU{2)r x U{1)b-l 
— > 5C/(3) X SU{2)l X U{1)y 
— > 5C/(3) X [/(I) 

EM 

(2, 

The U(2) family symmetry is broken in two steps H 
and the mass hierarchy is produced using the Froggatt- 
Nielsen mechanism [|| : 



Ui2) ^ Uil) nothing 



(3) 



and M is the UV-cutoff of the effective theory above 
which the family symmetry is exact, and eM and e M 
are the VEV's accompanying the flavon fields given by 
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The various aspects of VEV's of Higgs and flavon fields 
are given in CM. 

The superpotential for our model is 



a 



(5) 



Analytical Analysis of Mass Textures: In our 

convention, the Yukawa couplings F^, (i = u,d,e,ui^]i) 
are defined in the way that the left-handed fields are on 
the right, and the charge conjugate fields are on the left. 
In the quark sector, they read 



Woirac = V'sV'sTl + —Ip^ll^a {T2(j){l) + 33(/)(2)) 
+ J^'4'a'4'b {Ti + C) S'(2) + -^Ipa^JbTsS^i) 
W,^^ = iPs^sCl + ^7^3^a$C2 + jj^a^b^C2 (6) 

where Wfiavon is the superpotential involving only flavon 
fields which give rise to their VEV's given in Eq. (34) of 
CM Q. The mass matrices then can be read from the 
superpotential to be 

(lOj)e' 

Mu,.,n = I (10+) £ (l03+)e 
(10+).' (IO3+). (10+) 

rae' 

r4e e I Mu (7) 
r2e' e 1 



M, 



d.e 



(105-y 

(10.^) e' (I,-3)(f26" 



e' 
e' (I,-3)pe I Md 
I 




(8) 



where Mu = (IO+), Md = (lO^), ^2 = (l0+)/(l0+), 
r4 = (10|) / (lO;,^) and p = (j^^) I (lOr)- The right- 
handed neutrino mass matrix is 

/ 













(T26;°) 


Si] 











)52 


(T26;°) 


S3 




(1262° 




(1262°; 




(T26;° 


) ) 





82 

5i 63 




Mr 








(9) 



with Mji = (^126i'^. (This is one of the five sets of sym- 
metric texture combinations (labeled set (v)) proposed 
by Ramond, Roberts and Ross 0.) Here, the super- 
scripts +/ — /O refer to the sign of the hypercharge. It 
is to be noted that there is a factor of —3 difference be- 
tween the (22) elements of mass matrices Md and M^. 
This is due to the CG coefficients associated with 126; as 
a consequence, we obtain the phenomenologically viable 
Georgi- Jarlskog relations j|] . 

In CM, VEV's were taken to be real. In general, all 
VEV's are complex leading to spontaneous CP violation 
which is the subject matter of this paper. 



C-mass — —YuU rQlHu — YdDnQLHd + h.c. 
And the charged current interaction is given by 

Here we have written the Lagrangian in the weak basis. 
The Yukawa matrices are diagonalized by the bi-unitary 
transformations 



diag{yd,ys,yb) 



(10) 



where Vr and Vl are the right-handed and left-handed 
rotations respectively, and all the eigenvalues y^'s are 
real and non-negative. To extract the left-handed (right- 
handed) rotation matrices, we need to consider the di- 
agonalization of the hermitian quantity Y^Y {YY^}. 
The Cabbibo-Kobayashi-Maskawa (CKM) matrix is then 
given by 



Vud ^us ^ub 
VcKM^Vu^Vl^ ( V,d Vcs V,b 

Vtd Vts Vtb 



(11) 



The unitary matrix Vckm has in general 6 phases. By 
phase redefinition of the quarks, one can remove 5 of 
them. The remaining one phase is one of the sources for 
CP violation. A parameterization independent measure 
for the CP violation is the Jarlskog invariant ||] , defined 
as J^p = Im{ViiV^2^2i^22} ■ And the three angles of 
the CKM unitarity triangle are defined as 



1 = Arg{- 



^) 



(12) 



In the lepton sector, the charged lepton Yukawa matrix 
is diagonalized by 



ydraa ^ v,,Y,Vl = diag{y,,y^,yr) 



(13) 



And the effective light left-handed Majorana neutrino 
mass matrix (obtained after using the seesaw mechanism) 
is diagonalized by an unitary matrix V^^^ : 



lyfdiag _ 



V-^LLMuLLyuLL = dittg {m^^ , m^^ , m^^) (14) 



where the eigenvalues ye,n,T and mi,-^^^^^^^ are real and 
non-negative. The leptonic mixing matrix, the Maki- 
Nakagawa-Sakata (MNS) matrix, can be parameterized 

by 



3 



C12C13 S12C13 Si3 

UmNS = Ve^Vj^^ ^ \ -S12C23 - Ci2S23Sl3e**' -C12C23 - Si2S235l3e"'' 3230136"^' | ( ) (15) 

S12S23 - Ci2C23Sl3e*'^' -C12S23 - 5122233136"^' 2230136"^' 



Note that the Majorana condition, 



(16) 



where C is the charge conjugate operator, forbids the 
rephasing of the Majorana fields. Therefore, we can only 
remove 3 of the 6 phases present in the unitary matrix 
Umns by redefing the charged lepton fields. Note that 
Umns is the product of an unitary matrix, analogous 
to the CKM matrix which has one phase (the so-called 
universal phase). Si, and a diagonal phase matrix which 
contains two phases (the so-called Majorana phases), 0:21 
and a3i. The leptonic analog of the Jarlskog invariant, 
which measures the CP violation due to the universal 
phase, is given by 

fcp EE Im{U^2Ue3U;^U:^} (17) 

For the Majorana phases, the rephasing invariant CP vi- 
olation measures are ll(| 



Si EE Im{UelU:-i}, S2 EE Im{Ue2U:3} 



(18) 



From S'l and S'2, one can then determine the Majorana 
phases 

cosa3i = l-2 |^^^|f|V^^|, ^^g^ 
cos(a3i - Q21) = 1 - \ u,2nu,3\-^ 

The Lagrangian is invariant under the following weak- 
basis phase transformations 

Ur^u'j,^ KIUr, Dr -^d'j,^ K^Dr (20) 

Yu^Y;^^ K^YuKl, Yd-^Y'^^ K^YuKl (21) 

where K'l, K^, and Kf^ are diagonal phase matrices. We 
are interested in complex symmetric textures resulting 
from the 5*0(10) relations: 
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Y, 
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with a ~ 6 < c < 1 and e < / < 1. The weak- 
basis phase transformations mentioned above enable us 
to reduce the number of phases to two, 6 = (7^ — 27c — 7) 
and ^ = (7e — 7/ -l- 7c — 7a)- Then we have 



Y,, 




Y, 



ee-^« 
ee'« (1,-3)/ 1 h (22) 
1 



We then diagonalize the mass matrices analytically. In 
the leading order, the mass eigenvalues are given by 



C2| 



rric — |(— fee* +c )|du„ 
mt ~ dvu 

e'fh 



rrid 



e^ + f 

o2 



(2e^/ + /> 
m., ~ — Vd 



P 



nib = hvd 

3e^fh 



me 



9f 



;Vd 



6e^ fh + 27 fh 

ttIt- = hvd 
The CKM matrix elements are given by 



- + -e 
/ c 

Vcd ^ e'"' - - 

c / 



C 

e 



Vu ^ jc 
Vts =^ -c 



(23) 



(24) 



4 



The phase Sq is given by 

6q ~ tan-i(--sin(6l' -0) 



c 

e a 



tan-^(™--sin(^^ -O) + (0 -0 (25) 
/ c 



where 



6*' ~ tan~^(6sin6'- 2- 



■ sm 



6 

2— sin 



0) (26) 



The form of the CKM matrix in Eq. ( p^ is also the one 
favored in [Tl|| . 

In the lepton sector, similar diagonalization is carried 
out, and the charged lepton diagonalization matrix is 



Ve ^Vd{b^ -36) 



(27) 



rric = 677 MeV{677tlt) 

rat = 182 Gey(181^13) 

me = 0.485 Mey(0.486847) 

= 103 Mey(102.75) 

mr = 1.744 06^(1.7467) 

IK^I = 0.222(0.219- 0.224) 

iKbl = 0.0039(0.002- 0.005) 

\Vch\ = 0.036(0.036- 0.046) 

where the values extrapolated from experimental data 
are given inside the parentheses ||l^. These values cor- 
respond to the following set of input parameters at the 
GUT scale, Mqut = 1.03 x lO^^ GeV: 



a = 0.00246, b = 3.50 x 10" 



c = 0.0320 



We obtain bimaximal mixing in the neutrino sector and 
Am2j„ > Am| by choosing § 



d = 0.650, 61 = 0.110 




(28) 



In the basis where the Yukawa matrices take the forms 
of Eq. (IH), this implies that the right-handed neutrino 
mass matrix, Eq. (^, has the following elements: 



-^1 

63 



2a — 2ac ^ 



cH 



2a - 2ac + cH 
a{c - 6e^^) -I- bete" 
2a - 2ac + cH 



(29) 



The three eigenvalues of M^^^ are 



^V2 ^ Mr 



nil. 



(2 + -)- 



(30) 



Note that since Ve^ is approximately an identity matrix, 
Umns ~ Kxi,- Consequently, with the present experi- 
mental status, it is not possible to make unique predic- 
tions for the leptonic CP violating phases from the fitting 
of the absolute values of the MNS matrix elements, as 
was done in the quark sector. We therefore assume that 
the light neutrino mass matrix to be real. And the lep- 
tonic CP violation will solely due to the phase present in 

Ye. 

RGE Analysis: We use the following inputs at Mz ~ 
91.187 GeV: 



nil, 



2.32 MeV{2.33_ 



-0A2\ 
-0.45) 



e = 4.03 X 10~^ / = 0.0195 



h = 0.0686, e = -0.720 



91= 92 = .93 = 0.746 



(31) 



(32) 



the one-loop renormalization group equations for the 
MSSM spectrum with three right-handed neutrinos |l^ 
are solved numerically down to the effective right-handed 
neutrino mass scale, Mr. At AIr, the seesaw mechanism 
is implemented. We then solve the two-loop RGE's for 
the MSSM spectrum ||ll down to the SUSY breaking 
scale, taken to be mt{n^} = 176 A GeV, and then the SM 
RGE's from nit{mt) to the weak scale, Mz- We assume 



that tan/3 



10, with 



(246/V2 GeV)^. 



At the weak scale Mz, the predictions for 



9i 
4ir 



are 



ai= 0.01663, a2 = 0.03374, = 0.1242 

These values compare very well with the values ex- 
trapolated to Mz from the experimental data [p^ , 
(ai, 02,03) = (0.01696,0.03371,0.1214+0.0031). The 
predictions at the weak scale Adz for the charged fermion 
masses, the CKM matrix elements and the CP violation 
measures, are summarized in TABLE || (after taking into 
account the SUSY threshold correction ||l^, = —0.15) 
along with the experimental values extrapolated to Mz 

The current result from the atmospheric neutrino ex- 
periment is Jlql 



Am^3 = (1.3 - 8) X 10"3ey2(best fit : 3.3 



X lO"-*) 



sin2 26' = 1 



In the solar neutrino sector, the current best fit from the 
solar neutrino experiment data for each of the solution is 
given by 0: 
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TABLE I: The predictions for the charged fermion masses, 
the CKM matrix elements and the CP violation measures. 





experimental results 


predictions at Mz 




extrapolated to Mz 




rris 
i^d 


17~25 


25 


TTis 


93At\llMeV 


85.66Mey 


Til u 

• • I'D 


3 00'*'0 llGeV 


3 147GeV 


\Vud\ 


0.9745 - 0.9757 


0.9751 


\ ' ca \ 


0.218 — 0.224 


0.2218 


\Vcs\ 


0.9736 - 0.9750 


0.9744 


|Vtd| 


0.004 - 0.014 


0.005358 


\Vts\ 


0.034 - 0.046 


0.03611 


\Vtb\ 


0.9989 - 0.9993 


0.9993 


J'' 


(2.7ltl.l2) X 10"^ 


1.748 X 10"^ 


sin 2a 


-0.95 - 0.33 


-0.8913 


sin 2/3 


0.59t0.14t0.05 (BaBar) 
0.99t0.14t0.06(Belle) 


0.7416 
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34° - 82° 


34.55° (0.6030rad) 



LMA : (large angle MSW with larger mass squared) 
Am2 = 3.2 X 
tan2 e = 0.33, sin^ 261 = 0.75 
LOW : (large angle MSW with smaller mass squared) 
Am2 1 X IQ-'^eV'^ 
tan^ e = 0.67, sin^ 261 = 0.96 
SMA : (small angle MSW) 

Am2 = 5.0 X 10-^6^2 
tan2 e = 0.0006, sin^ 261 = 0.0024 
QVO : (quasi-vacuum oscillation) 

Am^ = 8.6 X 10~i°ey2 
tan2 6'=1.5, sin^ 261 = 0.96 

At present, none of these solutions can be ruled out p8| . 
The current bound on \ Uev3 \ from CHOOZ experiment is 

El 

|f/e.3l <0.16 

Our model favors both the LOW and QVO solution. The 
LOW solution is obtained with the following set of pa- 
rameters: 

Mr = 2.012 X lO^^GeV, t = 0.088 (33) 

which give rise to, using Eq. (|29|), 

5i = 0.001247 

82 - 2.221 X lO-^e* 

^3 = 0.01648e-'("°o"") 

This gives rise to three light neutrino masses 

m^, = 0.001711 ey 



= 0.001762 eF 
TO^3 = 0.05438 ey 

and the squared-mass differences are 

Am^,„ = 2.954 x lO'^ eV^ 
Am^ = 1.769 X 10"^ eV"^ (34) 

And the MNS matrix is given by 

/ 0.6743 0.7346 0.07497 \ 
\Umns\ = 0.5427 0.4322 0.7202 

\ 0.5008 0.5230 0.6897 / 

This translates into 

sin^ 29atm = 4|(7^.3l'(l - If^^-,,!") - 0.9986 
sin^ 200 = ^Uei,^^{l - iC/ei^^n = 0.9937 

and the MNS matrix element If/ei/sl is predicted to be 
0.07497, in very good agreement with the experimental 
result. The leptonic analog of the Jarlskog invariant is 
predicted to be 

fcp = Im{UiiU^2U2iU22} = -0.008147 

The matrix element for the neutrinoless double beta 
if^Poi') decay, | < m > |, is given in terms of the rephas- 
ing invariant quantities by 

|<m>|2 = ml\Uei\^ + ml\Ue2\^ + mlpesl^ 
+ 2mini2\Uelf\Ue2f cosa2i 
+ 2mim3|L/eiP|J7e3pCOSQ;3i 

-I- 2m2m3|L/e2nJ7e3pcos(Q;3i - a2l) 

(35) 

The current upper bound for | < m > | from the experi- 
ment is 0.2 eV |20|. The two Majorana phases (asi, a2i) 
are (0.9708, —1.326); they give rise to a highly suppressed 
I < m > I = 1.359 X 10-3 eV[^. The three heavy neu- 
trino masses are given by 

Ml = 9.412 X lO^Gey 
M2 = 1.486 X lO^Gey 
M3 = 2.013 X lO^^GeV 

The QVO solution is obtained with: 

Mr = 1.217 X t = 0.0142 (36) 

which give rise to, 

5i = 0.001267 

82 = 3.641 X lO-^e' ("-^^OD 

(53 = 0.01502e-'("°i°^4) 

The three light neutrino masses are predicted to be 

m^i = 2.856 X lO"'^ eV 
= 2.870 X 10"-^ eV 
m^, = 5.587 X 10"^ eV 
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and the squared-mass differences are 



Ar 



Ami 



= 3.122 X 10"^ 

= 7.584 X 10"^" eV^ 



And the MNS matrix is given by 

/ 0.6794 0.7318 0.05310 
\Umns\ = 0.5302 0.4513 0.7178 
(^ 0.5072 0.5107 0.6942 

This translates into 

sin^ 20at,n = 0.9991, sin^ 29^ ^ 0.9950 

and the MNS matrix element lUei^^l is predicted to be 
0.05310, in very good agreement with the experimental 
result. The leptonic analog of the Jarlskog invariant is 
predicted to be 

f^p = -0.008110 

The two Majorana phases (0131, a2i) are (1.385, —0.4976) 
and the matrix element for /3/3oi^ is predicted to be 
I < m > I = 3.07 X 10""* eV. The three heavy right- 
handed Majorana neutrino masses are 

Ml = 3.699 X lO^GeV 
M2 = 2.341 X 10^°GeV 
M3 = 1.218 X 10^'*GeF 

Our model can also produce the LMA solution by prop- 
erly choosing the values for Mn and t; however, the pre- 
diction for C/ei/3 violates the experimental bound, leading 
to the elimination of the LMA solution in our model. 

A few words concerning baryonic asymmetry are in 
order. Even though the sphaleron effects destroy bary- 
onic asymmetry, it could be produced as an asymmetry 
in the generation of {B — L) at a high scale because of 
lepton number violation due to the decay of heavy right- 
handed Majorana neutrinos which in turn is con- 
verted into baryonic asymmetry due to sphalcrons. But 



in our model this mechanism produces baryonic asym- 
metry of 0(10""'^'^) which is too small to account for the 
observed value of (1.7 — 8.3) x 10^^^ [Q, reasons be- 
ing that the mass of the lightest right-handed Majorana 
neutrino is too small and the 1 — 3 family mixing of right- 
handed neutrinos is too large, leading, in essence, to the 
violation of the out-of-equilibrium condition required by 
Sakharov [Q. So a mechanism other than leptogenesis 
is required to explain baryonic asymmetry. 

Summary: We have generalized our recently con- 
structed model based on SUSY 50(10) combined with 
[7 (2) family symmetry to include phases in the mass ma- 
trices leading to CP violation. In contrast with the com- 
monly used effective operator approach, 126-dimensional 
Higgs fields are utilized to construct the Yukawa sector. 
i?-parity symmetry is thus preserved at low energies. The 
symmetric mass textures arising from the left-right sym- 
metry breaking chain of 5*0(10) give rise to very good 
predictions for quark and lepton masses and mixings. 
The prediction for sin 2/3 agrees with the average of cur- 
rent bounds from BaBar and Belle. In the neutrino sec- 
tor, our predictions are in good agreement with results 
from atmospheric neutrino experiments. Our model al- 
lows both the LOW and QVO solutions to the solar neu- 
trino anomaly, and the matrix element for pPov decay is 
highly suppressed. The leptonic analog of the Jarlskog 
invariant, Jqp, is predicted to be of 0(10^^). It is in- 
teresting to note that, in the Yukawa sector, the model 
predicts (12-1-3) masses, (6 + 3) mixing angles and (4-f 3) 
phases (- the additional 3's in the parentheses refer to the 
right-handed neutrino sector) in terms of nine parame- 
ters given in Eq. (^ij) and t and Mr, a total of eleven 
parameters. 
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